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Abstract 

The form factors of the weak currents, which appear in the semileptonic decays of 
the heavy pseudoscalar mesons are calculated within the quark confinement model 
by taking into account, for the first time, the structure of heavy meson vertex and 
the finite quark mass contribution in the heavy quark propagators. The results are 
in quite good agreement with the experimental data. 



1 Introduction 



The study of semileptonic decays of heavy pseudoscalar mesons can be used to 
determine the elements of the Cabibbo-Kobayashi-Maskawa (CKM) matrix. 
The decay D K{K*)liy is related to \Vcs\, B D{D*)lu and B 7r(p)/z/ 
are proportional to IK^P and | V^bp, respectively. In the charm sector, however, 
the CKM elements can be determined independently of the D semileptonic 
decay rate using unitarity of the CKM matrix and the smallness of Vcb and 
Vub [1^3]. Thus, the theoretical predictions for the form factors and their q^- 
dependence can be tested. 

The study of heavy-to-heavy transitions in decays of -B — D{D*)lv is con- 
siderably simplified by the spin- flavor symmetry [4]. In the limit of infinite 
quark mass, in fact, the quark mass and spin decouple from the dynamics of 
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the decay, leading to numerous symmetry relations among form factors which 
can all be related to a single universal form factor, the Isgur-Wise function. 
At zero recoil, this function is known to be normalized to unity, which allows 
one to determine \Vcb\ from the measured B — > D*li> spectrum in the small 
region near the zero recoil point. The theoretical symmetry corrections are of 
1/mQ order due to the Luke's theorem [5]. 

The determination of \Vuh\ from analysis of S ^ 'k{p)Iv decays is one of the 
most important and challenging measurements in B-physics since the rate for 
these decays is expected to be only about 1% of the inclusive semileptonic 
decay rate. The exclusive calculations for B — > Xylv are more difficult than 
those for B — > XJv, because the range of recoil velocities available to the 
light final-state mesons is much larger than for the charm mesons. One there- 
fore expects a much larger variation in the form factors, which are still poor 
known, that enter into the decay rate. As a result, measurements of are 
currently quite model dependent, and there is substantial variation among 
values obtained using different models [2,3]. 

The main goal of the present paper is to describe the hcavy-to-heavy and 
heavy-to-light transitions within the quark confinement model (QCM) [6], by 
taking into account for the first time the nonlocal heavy-light quark vertices. 

The QCM approach is based on modelling the confined light quarks with 
the assumption of local hadron-quark coupling. It successfully describes many 
static and non-static properties of light hadrons. The extension of this ap- 
proach to heavy quark physics has been done in [7] by assuming that the free 
Dirac propagators can be employed for charm and bottom quarks. It might be 
justified by the observation that heavy quarks weakly interact with vacuum 
background fields, and therefore they can be considered as free particles with 
large constituent masses. The scaling laws for leptonic decay constants and 
semileptonic form factors are reproduced in the heavy quark limit. In addition, 
the Isgur-Wise function has been calculated. However, the Isgur-Wise func- 
tion is larger than in other approaches and in the fitted experimental data. In 
[8,9] the infrared behavior of the heavy quark has been taken into account by 
modifying its conventional propagator in terms of a single parameter v and 
the heavy-to-light form factors have been calculated. In this paper we intro- 
duce the vertex function describing the distribution of constituents inside a 
heavy meson. Such distribution is related to the heavy-meson Bethe-Salpeter 
amplitude in the approach based on the Dyson- Schwinger equations [10] 
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2 Model 



The QCM approach [6] is based on the effective interaction Lagrangian for 
the transition of hadron into quarks: 

Cint{x) = gHH{x) jdXijdX2^H{x]Xi,X2)q{xi)VH\Hq{x2) ■ (1) 



Here, \h and Th are the Gell-Mann and Dirac matrices, respectively, which 
provide the flavor and spin numbers of mesons H. The function is related 
to the scalar part of Bethe-Salpeter amplitude. The local form ^h{x] Xi,X2) — 
8{x — [xx + X2)l'2)b{x\ — X2) has been used in the QCM [6]. 

The coupling constants qh defined by what is usually called the compositeness 
condition proposed in [11] and extensively used in [6], is given by 

Zh^i- §n^«) = (2) 

where 11'^ is the derivative of the meson mass operator. 

In the QCM-approach the light quark propagators are given by an entire (non- 
pole) function to ensure the quark confinement: 



niq— J A//— 



f da,, „ 1 r . p\ ^ P% 



(3) 



with the functions a and b defined by 

„(-.) = fJp^ 6(-,) = [4^. (4) 

J U'^ — z J — z 



Moreover, to conserve the local properties of Feynman diagrams like the Ward 
identities, one has the following prescription for the modification of a line with 
n-light quarks within the Feynman diagram [9]: 

n;^r.^/..„nx^r.. (5) 
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Zl - Z2 

b{z,) - b{z2) 



It is useful to introduce the notation 

jjj ~r Z jJj ~r Z 

j daf,as{z) = a{z) J da^av{z) = b{z) 
J dau,as{zi)av{z2) = J dau,av{zi)as{z2) = ^^^^^ "'^^'^^ 
/ daf,av{zi)av{z2) - 

J Zl — Z2 

f f \ f \ zib{zi) - Z2b{z2) 
/ da^as{zi)as{z2) = 

J Z\ — Z2 

J Z\ — Z2 

where the confinement functions employed in [6] have the forms: 

a{u) — oo exp(— — a\u) b{u) = bo exp{—u^ + biu) . (6) 

The following values for the free parameters Oj, 6j, and A: 

ao = 60 = 2 ai = 1 61 = 0.4 A = 460 MeV; 

give a good description of the hadronic properties at low energies [6] . 

The hadron-quark coupling constants for light, pseudoscalar and vector, mesons 
and heavy pseudoscalar mesons are also determined from the compositeness 
condition [6] and written down 

27r I 2 „ / X „ a; }, ux.(l-u/2) 

9p^^\hr-. r, Rp(x) ^ Bo + - dub( — -) , ' -. (7) 







Note that from now on all masses and momenta in the structural integrals are 
given in units of A. 

The heavy quark propagator is given by 



3 Form factors 



We consider the leptonic H[p) — > li>, semileptonic heavy-to-heavy B{p) — > 
D{p')lv and semileptonic heavy-to-hght H{p) P{p')lv decays, where H{p) 
represents a B (or D) meson with momentum p {p^ — m\j) and P{p') can 
be a TT or X meson with momentum p' {jP — m^). The invariant amphtudes 
describing the decays are: 



A{H{p) eu) = ^VQ,{eO,u)MM (9) 
A{B{p) ^ D{p')ei^) = ^VUeO,iy)M^^{p,p') (10) 
A{H{p) P{p')eu) = ^VQ,{eO,u)M^p{p,p'), (11) 

where Gp is the Fermi weak-decay constant, Vgq is the appropriate element 
of the Cabibbo-Kobayashi-Maskawa matrix (g denotes a light quark and Q a 
heavy quark) and the matrix elements of the hadronic currents are: 



MM 



47r2 
xtr 



(12) 



xtr 



(13) 



xtr 



0'^'5q(^+:^)7V(^)7V(^+^) 



^fni")ip+p'r+f-'^iQ')ip-p') 



(14) 
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Prom the compositeness condition (in Eq. (2)), the expression for the propa- 
gators, in Eqs. (3), (4) and (8), and the method outhned in [10], we obtain for 
the heavy decay constants and heavy to heavy form factors 



9h 



47^2 



\ 3Jf\mH,mH) 



47r2 



9h J2{mH) 



fBD 
J± 



47r2 



(15) 



with 



/u 
dUj^-^—:^z'(f)H{z) 




X 



MQa{z) + -b{z) (2z + u(mjj + + z)) 



1 CXJ 

/r'(mB,m«) = i/rfx/ 



^ 1 oo 

- y" rf^y" du-—^—^(j)B{z^)(j)D{z^) {a{z:c) {Mh + Mc) 

0^ ^ 

+b(z^) [u [MbM^ + ml(l-x)+ xm\ + z^) + 2z^] } 



1 oo 



4 \mB,mD) = ^ J dx J du—^—-4>B{z^)(pD{zx) 



0^ ^ 

{a{z^) [Me -Mb + 2u (Me - x{Mb + M^))] 

+b(z^)u [(1 - 2x)(z^ - MbM^) + ml{l - x) - xm|] } 



where the variables are given by 



z = uM?, 



u 



1+U ^ 



z'^M% 



(l + « 



Za- — u <x 



m 



B 



1 + u 



+ {l-x) 



Mi 



m 



D 



1 + u 



u 



l + u 



x{l — x)(f 



For the heavy to hght form factors, instead, the analytical expressions are 
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47r2, 



2 oo oo ^ 



^ ^ _1 V ; 



where the functions Gi{zi, z-i) and G2(-2i, -^2) can be written as 



G2^Mq{1 + u)Fsv(zi, Z2) + {zi(l + u)+ umjj) Fyy(zi, Z2) 

+2t^ (Fss'izi, Z2) + ZiFyy>{zi, Z2)) 



and 



= + uMn 



Z2^X2 + iy2 = 



um 



H 



1 + u 



2 , ..2 m 

r + uMq - - 



m. 



+ u 1 + u 



r 



+ i 



2r'~fmp 



The functions Fu appearing in Gi and G2 are defined as: 
Fss{zi,Z2) = j danas{zi)as{z2), Fyy>{zi, Z2) = j da,^ay{zi)ay{z2), etc. 

Before closing this section, we discuss the behaviour of the heavy-to-heavy 
form factors in the hmit of Mi,, Mc — > 00. We shall show that our model 
reproduces, in this limit, all the scaling laws predicted by the Heavy Quark 
Effective Theory at leading order. 



In particular, in the heavy quark limit [mjj — {Mq + EY and Mq — > 00] one 
finds 



3^ J_ 
47r2 2Mq 



Mq 27r V Ihh 



00 

Ihh = j du(j)]j{z){a{z) + ^/uh{z)}, 



du{y/u - E)(l)H{z){a{z) + ^^/ub{z)}, 



2^MqMc 



Q' 
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where the Isgur-Wise function, ^{w), is given by 

C{w) = / 777 / du(t)]j{zw) a{zw) + Ju/Wh{zw) 

Ifjff J W J L 







(16) 



with 



W = 1 + 2t{1-t){w -1) zw = u-2EJu/W 



z = u 



It is readily seen that the upper bound for the Isgur-Wise function is obtained 
for £^ = 0, namely 



1 nn[w + Vw^ - 1] 2R 
1 + r\ Vw2-1 ^ 1 + w 



(17) 



where 



R 



J du 4>h{u) y/u b{u) 



oo 

/ du (f)H{u) a{u) 



As a consequence of Eq. (17) the slope parameter has the lower bound 



= -cm = J 



21 + R. 



1 

> - 
- 3 



In the heavy quark limit {p^ = {Mq + E)"^, [p'Y = and Mq — >■ oo) one finds 
for the heavy-to-light form factors that 

1 



Mr 



.(18) 



Here 



Gi = Fss{zi, Z2) + ziFvv{zi, Z2) 

G2 = Fsv{zu h) + Ty/uFvvih, Z2) 
with the F/z's defined before, zi — u — 2E^Ju, Z2 — Z1 + 2Xt^Ju, and 



Mr 



Ml 



At the end point — {X — 0) one can reproduce the well-known relations 
among form factors in the heavy quark limit 



4 Results and discussion 

The expressions obtained in the previous section for the form factors and decay 
constants are valid for any kind of vertex function /c^). Here, we choose a 
Gaussian form 0(— /c^) = exp{A;^/A^} in Minkowski space. The magnitude of 
characterizes the size of the BS-amplitude and is an adjustable parameter 
in our approach. Thus, we have four adjustable parameters: Ad and Kb plus 
the two heavy quark masses, or binding energies Ed — uid — Mc and Eb — 
niB — Mb. The first two are fixed in such a way the form factors /^^{q^) 
and f^^{q^) are increasing functions of q"^; we choose Ad = 0.56 GeV and 
Ab = 0.67 GeV. The other parameters are fixed by the least-squares fit to the 
observables measured experimentally or taken from a lattice simulation (see 
asterisks in Table 1). 

The best fit is achieved for Ejj ~ Eb, thus we choose to fix E^ = Eb in such a 
way we have only two free parameters. The best values are E^ = Eb = 0.554 
GeV and Vcb = 0.043 which is close to the world-accepted value [1]. The 
resulting values for the heavy to light form factors at = are larger those 
predicted by other approaches. It should be stressed that these values are 
practically fixed by the assumption that the form factor should be increasing 
functions of q'^. Moreover, there is a strong correlation between /:^'^^(0) and 
the decay constant fn, i-e. smaller values for form factors corresponds to small 
values for decay constants. The situation changes if no assumptions are done 
on the q'^ behaviour of the form factors. 

We plot the the g'^-behaviour of the resulting form factors on Fig.l. For com- 





9 



parison, the vector dominance, pole model is shown: 

QQ 

with rriY being the mass of the lightest gg'-vector meson. We use mo* — 2.11 
GeV for c ^ s, uib* — 5.325 GeV for b ^ u, uib* ~ — 6.4 GeV [12] for 
b ^ c transitions. The values of (0) are taken from the Table 1. Also we 
calculate the branching ratios of semileptonic decays by using widely accepted 
values of the CKM matrix elements [1]. 

A few comments should be done concerning the comparison of our results with 
the results of paper [13] where the weak decays of pseudoscalar mesons have 
been described within the relativistic constituent quark model with free quark 
propagators. Since there is no confinement in that model the binding energies 
have been found to be relatively small: Ed = 0.20 GeV and Eb = 0.22 GeV. 
Such values provide the absence of imaginary parts in the physical amplitudes 
describing the decays of the low-lying pseudoscalar mesons. However, the ex- 
cited states like vector mesons cannot be considered in a self-consistent man- 
ner. The Quark Confinement Model allows us to give the unified description of 
physical observables without quark thresholds in the physical amplitudes and 
with a minimum set of parameters: the only parameter A = 0.460 GeV, the 
size of confinement region, for light quark sector and four extra parameters 
(As^ij-the sizes of Bethe-Salpeter amplitudes, and i^^^D-the binding energies) 
for heavy quark sector. As a result, the accuracy of desciption is less than in 
[13] while, the region of application is considerably wider. 
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Table 1 

Prediction for leptonic decay constants (in GeV), form factors and ratios. The 

(niB - mof)- 
Obs. Calc. Obs. 



"Obs." are extracted from Refs. [1,14-19] {q\j 



Calc. 



* fo O.I9IT28 0.165 

* /+^(0) 0.74 ± 0.03 0.77 

* I Kb I /f^ (9m) (5.09 ± 0.81) 10-2 5.110-2 

* /f ''(O) 0.27 ± 0.11 0.55 



0.172^^^ 



0.135 



* /b 

Br(L> Klv) (6.8 ± 0.8) • 10"^ g g . 

Br(S Dlu) (2.00 ± 0.25) • IQ-^ 3.5 • 10-^ 

Br(S ^ ttIv) (1.8 ± 0.6) • lO"'' 3.3 • 10""^ 



contract number 99-02-17731-a. 
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V.t=0.0032 

Br(B - 7T I l/)=3.3 10"' 

Br(B - 71 I !/)„,=( 1.8 +/- 0,6) 10"' 




q^(GeV^) 



Fig. 1. The semileptonic D ^ K, B D and B ^ ir form factors with, for 
comparison, a vector dominance, monopole model Eq. (19) and a lattice simulation 
[20]. Our results: solid lines. Monopole: dotted lines. Lattice: data points. 
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